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00 . Abstract 

cn 

^^ ' We study the production of gravitinos during the preheating era af- 

^P ■ ter inflation by means of the non-perturbative Bogolyubov technique. 

Considering only the helicity ±3/2 states, the problem is reduced to 
the simpler Dirac fermions case. We calculate the production in a par- 



o ■ ticular supergravity model in an expanding universe and obtain the 

p. 



f^ . spectrum and number density. Finally we compare the results with 

D I the nucleosynthesis bounds and extract some consequences. 



1 Introduction 

The presence of weakly interacting massive particles can have important 
consequences in cosmology [|I|]. In particular, when such particles are light, 
with life times longer than the universe age, they can act as dark matter or 
even dominate the present energy density and overdose the universe. On 
the other hand, if they are unstable and decay during the nuclesynthesis era, 
they could destroy the nuclei created in this period and spoil the successful 
predictions of the standard big bang model. 
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One of such particles is the gravitino, i.e, the spin 3/2 superpartner of 
the graviton in supergravity theories, whose couphngs to the rest of particles 
are suppressed by the Planck mass scale. Typically its decay rate is given by 
T^I2 — m^,2/Mp which implies that gravitinos lighter than 7713/2 < lOOMeV 
will live longer than the universe age. Due to their weak couplings, graviti- 
nos freeze-out very early when they are still relativistic and therefore their 
primordial abundance can be estimated as 773/2/5 — 10~^ [Q. This allows us 
to obtain the well-known bound on the mass of stable gravitinos, i.e, in order 
not to overdose the universe we need 7773/2 < IkeV 0. However if they are 
unstable, that primordial abundance would give rise to an enormous amount 
of entropy that will conflict with the standard cosmology. A possible way out 
of this gravitino problem is the existence of a period of inflation that dilutes 
any primordial density 0. Unfortunately the problem can be recreated if af- 
ter inflation, gravitinos are produced by some mechanism. In fact, this could 
be the case if during the period of infiaton oscillations at the end of inflation, 
the reheating temperature was sufficiently high. A successful nucleosynthe- 
sis period requires (we give some conservative bounds |^): 773/2/s < 10~^^ 
for a gravitino mass 7773/2 — lOOGel^, 773/2/s < 10~^^ for 7773/2 — ITeV and 
'^3/2/s < 10^^^ for 7773/2 — lOTeV. The production of gravitinos during re- 
heating depends on the temperature Tr as 0: 773/2/s ~ IO'^'^Tr/ {10^ GeV) , 
which implies that for a typical mass 7773/2 — ITeV, the reheating temper- 
ature has to be Tr < lO^GeV. Another constraint appears in supergravity 
models, where the gravitino mass is determined by the scale of supersymme- 
try breaking. In order to solve the hierarchy problem, it is then suggested 
that 7773/2 < ITeV 0. 

As we have just mentioned, during the reheating period gravitinos can be 
created by the perturbative decay of other particles produced from the infia- 
ton oscillations. However, in recent years, the standard picture of reheating 
has changed dramatically as a consequence of some works in which it was 
realized that during the first infiaton oscillations, reheating can not be stud- 
ied by the standard perturbative techniques 0, |^. This so called preheating 
period can give rise to an explosive production of bosons due to the phe- 
nomenon of parametric resonance. In this period, the energy of the coherent 
oscillations of the infiaton field is very efficiently converted into particles. In 
the case of fermions, the limit imposed by Pauli exclusion principle avoids 
the explosive production, but however the results deviate from the pertur- 



bative expectations [^ [T0|, |IT] . This fact will be of the utmost importance 



in the case that gravitinos are directly coupled to the inflaton, since the new 
preheating period can give rise to an excess of production. 

In this work we present our results on the production of spin 3/2 particles 
from the non-perturbative Bogolyubov technique |jl2| and apply the method 



to the preheating period in a particular supergravity model. 

2 Quantization of spin 3/2 fields in external 
backgrounds 

Preheating is based on the phenomenon of creation of particles from classical 



backgrounds (see |T3[ and references therein). The spin 3/2 case is however 

special due to the consistency problems that avoid the quantization of such 

fields in the presence of scalar, electromagnetic or gravitational backgrounds 

1^. The only theory in which these problems seem to be absent is super- 



gravity, provided the background fields satisfy the corresponding equations of 
motion |T^. But, the complicated form of the Rarita-Schwinger (RS) equa- 
tion makes it very difficult to extract explicit results even in simple cases. 
Let us start by considering the massive RS equation in a curved space-time. 
We will include the coupling to a scalar field x by modifying the mass term, 



fwe will follow the notation in 16 



e^^P^^.^^Dpi,, + -(m3/2 - X)[7^71V'. = 0. (1) 

As usual in supergravity models we will consider Majorana spinors satisfying 
ip^ = CipJ^ with C = i'j'^'j^ the charge conjugation matrix. We have intro- 
duced the curvature of space-time by minimal coupling as done in supergrav- 
ity, i.e, Dpip„ = {dp + |r2^^Safe)^(j with fi^^ the spin-connection coefficients 
and T,ab = |[7a,7b]- The e'^^''" removes the Christoffel symbols contribu- 
tion in the covariant derivative. Following the analogy with the creation of 
gravitons from Einstein equations, we will consider only the linearized equa- 
tion in 1/M {Mp = SttM^) for supergravity 0, i.e, we will consider only 
the symmetric part of the spin-connection, ignoring the torsion contribution 
that is 0{M~'^). In flat space-time and when x = 0, the general solution 
of these equations can be expanded in four helicity / = s/2 + m modes 



(/ = ±3/2,±l/2): 

^p/(^) = e-ip- ^ J^^u{p, s)e,{p, m) (2) 

with Jsrn the Clebsch-Gordan coefficients whose values are: J_i_i = Jn = 1, 



J_ii = Ji_i = 1/v^ and J_io = Jio = \/ 2/3. u(p, s) are spinors with definite 
hehcity s = ±1 and e^(p, ?7i) with ?7i = ±1, are the three spin 1 polarization 
vectors. 

In order to study the preheating era, we need to consider the presence 
of the inflaton field and the expanding universe. With that purpose we take 
a scalar field only depending on time x(^) aiid a spatially fiat Friedmann- 
Robert son- Walker metric. 

Now, the expression in @ is not a solution of the equations of motion. 
However, we can look for general homogeneous solutions of the RS equation 
in the form: 

rj{x) = e'^^-'r\t) Y: Jsmuip, s)e,{p, m) (3) 

These fields satisfy the condition ''j^ip^ = and if we constraint ourselves to 
the helicity I = ±3/2 states, then the RS equation reduces to a Dirac form 



{lip- m3/2 + X)V^M = (4) 

and the condition D^ip^ = also holds. As far as these modes satisfy a Dirac- 
like equation, it suggests that all the difficulties in the gravitino quantization 
procedure would reduce to the helicity ±1/2 modes in this case. In fact 
the above ansatz (^ is not a solution for the helicity ±1/2 modes even for 
homogeneous backgrounds. 

Thus we see that the production of helicity ±3/2 gravitinos during pre- 
heating in an expanding universe can be studied in the same way as the 
production of Dirac fermions. With this purpose we have to reduce equation 
(^ to a second order equation. Let us first write the equation in conformal 
time defined as dt = a{ri)dri: 

(^ia~'-f^'d^-ms/2 + X + i^-^l'')i^i. = (5) 



where a = da/drj. We will take the following ansatz on the helicity / = ±3/2 
solutions: 



V'{x) = a-'l\r^)e^^^-'uf{r^) (6) 



with 



+ a{ri){m?,i2 - xiv)) fpi{v)u{p, s)e^{p, m) (7) 

and the normalization t/^^(0)[/^(0) = 20; where m-g^ = a(0)m3/2. It is possi- 
ble to check that this ansatz automatically satisfies 'y'^ipfj, = and D^ip^ = 0. 
An appropriate form for the spinor u{p,s) and polarization vectors en{p,m) 
can be obtained if we choose the Dirac representation for the gamma matrices 
and we take (without loss of generality) the 2-axis along the p direction. In 
this case uip, if = (1, 0, 0, 0), ^(p, -1)^ = (0, 1, 0, 0), e^(p, 1) = ^(0, 1, ^, 0) 
and e^(p, —1) = -75(0, 1, —i, 0). With this choice u{p, ±1) are eigenstates of 
7° with eigenvalues +1. Then equation (|^) reduces to the well-known form: 

ij^ + P"^ - i^ {a{v){m3/2 - Xiv) 

+ a\r]){my,-x{v)f)fAv) = (8) 

In order to quantize the modes we will expand an arbitrary solution with 
helicity / = ±3/2 as: 

<(^) = / (^^«"'^'(^) {^''-'U^^ivh,^ + e-''^-'U^^''{v)a^fl) (9) 

where the creation and annhilation operators satisfy the anticommutation 

relations {aj^,alfi'} = {2Ti)^2uj5iii5{p — p'). 

3 A supergravity inflation example 

We will consider a specific supergravity infiationary model (see ^Jl), in which 
the infiaton field is taken as the scalar component of a chiral superfield, 



and its potential is derived from the superpotential I = (A^/M)(0 — M)^. 
This is the simplest choice that satisfies the conditions that supersymmetry 
remains unbroken in the minimum of the potential and that the present 
cosmological constant is zero. CMB anisotropy fixes the infiationary scale 
around A = A/M ~ 10~^. For the sake of simplicity, we will consider the 
case in which the gravitino mass is much smaller than the effective mass of 
the inflaton in this model, 777,3/2 <C m^ ~ 10~^M and since the production 
will take place during a few infiaton oscillations, we will neglect the mass 
term in the equations. The scalar field potential is shown to be stable in the 
imaginary direction and therefore we will take for simplicity a real infiaton 
field. Along the real direction the potential can be written as: 

V{<P) = aV' ((2(0-1) +0(0-1)2)2-3(0-1)^) (10) 

where we are working in units M = 1. This potential has a minimum in 
= 1. The coupling of the infiaton field to gravitinos is given by the following 
term in the supergravity lagrangian 0: 

Ant = -^e^/2^^[7^7lV'. (11) 

e^/2 = AV/2(0-i)2 (12) 

where the Kahler potential has been chosen in such a way that the kinetic 
terms for the scalar fields are canonical G'($, $''') = ^^^ + log|/p. The 
infiaton and Friedmann equations can be written in conformal time as: 

+ 2^0+^V;^ = O (13) 

where the derivatives are with respect to the new time coordinate fj = aQX'^rj 
and the new scale factor is defined as b^fj) = 0(77) /ao with a^ = a(0). The 
solution of this equation shows that after the infiationary phase, the scalar 
field starts oscillating around the minimum of the potential with damped 
amplitude. Substituting in (|]) and taking x = e'^'^^ for this particular case, 
we obtain: 

_ J_(feO,.) ^ l_A fM) ^ (15) 




N, 



kI 




V 



Figure 1. -Number of helicity I = ±3/2 gravitinos (-/VkO as a function of time 
77 for K = 5 



with K = p/(aoA^). From this expression we see that when the scalar inter- 
action is switched off, even in the expanding background, there is no particle 
production. This is not surprising since in that case, the equations of motion 
(in the masless limit) are conformally invariant. Following ||^, |13] we can 
calculate the occupation number: 

N.i{f) = ^ (2/. + ^[t{f)Uf) - f:,{f)UT)] - ^&e^(^)/V./(T)P]) (16) 

In order for the particle number to be well defined, we will evaluate it 
when the interaction is vanishingly small, that is, for large values of T. Here 
Jki is a solution of equation (|15D with initial conditions /k/(0) = 1 and /k/(0) = 
—in which corresponds to a plane wave for fj < 0. In order to define the initial 
vacuum at 77 = 0, we have taken the inflaton to be at the minimum of the 
potential at that moment (0(0) = 1), which implies that the interaction 
term vanishes, i.e, e*^*-"^^^-''^ = 0. According to the definition of b we also 
have 6(0) = 1. We have chosen 0(0) = 1.8 in our numerical computations 
which corresponds to an initial amplitude of the inflaton oscillations around 
0.06M„. 
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Figure 2. -Number density of hclicity / = ±3/2 gravitinos (k^N/^i) against 

K. 

In Fig.l, the behaviour of N^i has been plotted as a function of time. We 
see how the production takes place in a few inflaton oscillations and that 
asymptotically the number of particles created is a well-defined quantity 
since the interaction vanishes. The results for the spectra in the expanding 
background can be found in Fig. 2. Notice that we have not considered the 
backreaction effect of the produced particles on the scalar field evolution. 
Because of the expansion of the universe, the production is significantly re- 
duced with respect to the fiat space case. As expected the resonance 
structure is affected by the expansion and the Pauli limit is not saturated, 
but the number of particles that are produced is not negligible. 

We expect modifications in the spectrum for different initial conditions, 
although the total number of particles created is not very sensitive to these 
changes. From Fig. 2, we can estimate a lower bound to the total number 
density of gravitinos with both helicities as: 



1 /■'^ ,. o , aiX 



n< 



3\6 



(^) = ^-^/ N,ip'dp=^l N^iK'dK (17) 



TT'^a^{r]) Jp^,„ ^ n'^a 



with pmin = ^TxHirf). Since today, (?7o,a(?7o) = 1) H « MaoX^, we get: 
npreivo) > a^lO'^^GeV^. To be compared with the number density of a ther- 
mal distribution of hehcity ±1/2 gravitinos as estimated in (the hehcity 
±3/2 could be even less dense): ^(r^o) — lO'^'^GeV^. The comparison de- 
pends on the dilution term ajj, i.e. on the scale factor at the end of inflation, 
and it shows that, for example, for a typical value |]l[ Oq — 10~^^, the vac- 
uum fluctuation production is suppressed with respect to the thermal distri- 
bution by a factor > 10~^°. The corresponding cosmological consequences 
have been studied in [If. Comparing with the entropy density today we get: 

If we compare these results with the nucleosynthesis bounds given in 
the introduction, we see that the preheating production is compatible with 
them for a mass of (unstable) gravitino 7713/2 > O{10TeV) for this particular 
model. This implies that gravitino masses lOOMeV < 7773/2 < lOTeV would 
be forbidden. On the other hand, if gravitinos are stable particles, then in 
order not to overdose the universe we get: 7713/2 < ITeV. These bounds 
should be taken with caution, as the results strongly depend on the model 
parameters. 

In the standard scenario of reheating for this model, the reheating tem- 
perature is very low Tr ~ lO^GeV which implies that the perturbative grav- 
itino production, that can take place either from 2 — > 2 proccesses involving 
gauge bosons and gauginos or from direct inflaton decay Jl^, is very small 



fireh/s ~ 0(10^^^). Comparing with the preheating results above, we see 
that the difference between the two mechanisms is apparent in these estima- 
tions. 
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